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Abstract 

Angular x-ray cross-correlation analysis (XCCA) is an approach to study the structure of dis- 
ordered systems using the results of coherent x-ray scattering experiments. Here, we present the 
results of simulations that validate our theoretical findings for XCCA obtained in a previous paper 
[M. Altarelli et al, Phys. Rev. B 82, 104207 (2010)]. We consider as a model two-dimensional (2D) 
disordered systems composed of non-interacting colloidal clusters with fivefold symmetry and with 
orientational and positional disorder. We simulate a coherent x-ray scattering in the far field from 
such disordered systems and perform the angular cross-correlation analysis of calculated diffraction 
data. The results of our simulations show the relation between the Fourier series representation 
of the cross-correlation functions (CCFs) and different types of correlations in disordered systems. 
The dependence of structural information extracted by XCCA on the density of disordered systems 
and the degree of orientational disorder of clusters is investigated. The statistical nature of the 
fluctuations of the CCFs in the model 'single-shot' experiments is demonstrated and the potential 
of extracting structural information from the analysis of CCFs averaged over a set of diffraction 
patterns is discussed. We also demonstrate the effect of partial coherence of x-rays on the results 
of XCCA. 

PACS numbers: 61.05.cp, 61.43.-j, 61.43.Dq, 61.43.Fs 
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I. INTRODUCTION 



The study of angular correlations in diffraction patterns has a long history with two 
main directions. On one hand, it has been used in small angle x-ray scattering experiments 
as a tool in the attempt to solve structures of molecules in solutions or, more generally, 
in non-periodic systems. It goes back to the work of KamP^, over 30 years ago; in much 
more recent times, the perspectives opened by ultra-short and ultra-bright pulses from free- 
electron laser^^, which provide an opportunity to acquire diffraction data in a time short 
compared to translational and rotational relaxation times, have rekindled interest in this 
approach!^. In these recent applications, the angular correlation is defined as an average 
over many diffraction patterns. 

An alternative purpose, which may seem more modest at first, but could be of paramount 
importance in the physics of disordered or partially ordered systems is the unveiling of hidden 
symmetries in a disordered collection of elements. This leads to the problem of understanding 
systems with the so-called medium range ordei'^^ as well as complicated dynamics and 
correlation between dynamical heterogeneity and medium-range order in a large class of 
glass-forming liquid^^^HE] such systems the relevant question will be, for example, can 
one recognize and identify an ra-fold symmetry axis of an individual molecular species from 
the diffraction patterns of a liquid composed of such molecules? Can bond angles be detected 
from the study of angular correlations of the diffracted intensity of an amorphous system? 
There are partly affirmative answers in the study of partially ordered quasi two-dimensional 
systems like liquid crystal^ in which hexatic bond order can be detected by the study 
of the angular dependence of the diffracted intensity. More recently Wochner et al^^^^ 
reported angular correlations with pronounced periodic character in a colloidal suspension 
of polymethylmethacrylate (PMMA) spheres, expected to form icosahedral clusters near the 
glass formation concentration. 

In our first paper'^, angular correlations in the diffraction patterns from a disordered col- 
lection of symmetric species were discussed theoretically in a general frame. The information 
content of angular correlation data in different experimental conditions (dilute versus dense 
systems, different statistics of orientational disorder, and so on) was derived. It transpires 
that generally the richest information is obtainable from dilute systems in which not too 
many molecules or clusters are found instantly in the illuminated volume. At the same 
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time for such systems, static fluctuations of the scattered intensity are relatively high, and 
do not give conflrmative information on the local structure of disordered systems. It also 
appears that partial orientation of the species, with lining up of the symmetry axes in one 
direction, enhances the effects considerably, in agreement with the results of other groups. 
Nonetheless, some aspects of the data by Wochner et alP^ were difficult to reconcile with 
our results, for example the pronounced fivefold symmetry of some of the observed angular 
correlation patterns. 

In this paper we apply the x-ray cross-correlation analysis (XCCA) to the study of the 
local structure (LS) of disordered systems. In general, the role of LSs can be played by 
clusters in the gas phase or a colloidal system, protein molecules, viruses, or complex bi- 
ological systems in solution. Here, we consider the particular case of model 2D systems 
composed of non-interacting clusters with identical shape. The main goal of this paper 
is to identify experimental conditions when symmetry and/or structure of clusters com- 
posing a disordered system can be correctly determined from the XCCA. By varying the 
density of clusters and characteristics of their orientational disorder, we analyze how these 
changes infiuence the angular cross-correlation functions (CCFs) and their Fourier spectra. 
We demonstrate a statistical behavior of CCFs calculated for different single realizations 
of a system and compare them with CCFs averaged over many realizations of the system 
(ensemble of diffraction patterns). Such analysis can be related to so-called 'single-shot' 
experiments, when the measurement time is much shorter than any typical structural relax- 
ation time of the system under investigation. Such type of experiments can be performed, 
for example, on fourth generation x-ray sources, such as free-electron laser^^ (FELs), that 
produce ultrashort pulses on the femtosecond scale. We also investigate the influence of 
partial coherence of the incident x-ray beams on the outcome of the XCCA. 

This paper is organized as follows. In the next section a brief theoretical description 
of the XCCA based on our previous work^*^ is given. In the third section a model of a 
coherent scattering experiment on a disordered system used in simulations is introduced. In 
the fourth section we study the CCF as a function of orientational disorder in the systems 
with a large number of clusters. Statistical fluctuations of the CCF in a model 'single-shot' 
experiment for completely disordered systems with a small number of clusters are analyzed 
in the flfth section. Results of calculations of the averaged CCFs for systems with different 
number of particles and degree of orientational disorder are presented in the sixth section. 
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In the seventh section we demonstrate the role of coherence properties of the incoming x- 
ray wavefields for the analysis of the local structure of disordered systems. The paper is 
completed by the conclusions and outlook section. Detailed derivations are summarized in 
the Appendix section. 



II. THEORETICAL BACKGROUND 



The intensity CCF Cq{A) as a function of the angular coordinate < A < 27r is defined 
as follow^ 

(/(g, ^)I{q, + A))^ - (J(g, y,)) J 

where I{q, if) is the intensity scattered at the momentum transfer vector q, is an angular 
coordinate around a diffraction ring of radius q, and {f{^))ip = 1/2-k J^^ f{(p)dip denotes the 
angular average around the ring. 

It is convenient to analyze the CCF Cq(A) using a Fourier series decompositiorP^, 

00 

C,(A) = 2^0"^ cos(nA), (2a) 

n=l 

C^ = - r C,{A)cos{nA), (2b) 
Jo 

where are the Fourier components. 

The Fourier components of the CCF Cg(A) for n 7^ are determined by the angular 
distribution of intensity /(g, (p) according to the following relatiorP^ 

C-=\nq)/I%q)\\ (3) 

where /"(g) = 1/ (27r) J^^ I(q, ip)e~^'^'^d(p are the components of the Fourier expansion of the 
intensity on the ring of radius q. According to its definition the zero Fourier component 
of the intensity distribution I^{q) coincides with the small angle x-ray scattering (SAXS) 
form-factor S{q). We also want to note here that normalization in Eq. ^ has to be done 
with caution, the g— values where /°(g) = should be excluded from the consideration. 

In the following we will consider a specific case of two-dimensional (2D) disordered sys- 
tems that on one hand simplifies XCCA analysis but at the same time reveals important 
general features common for such systems. It was shownP^, that in this case only even 



(n = 2/, / = 1, 2, 3, . . . ) Fourier components of the intensity I"'{q) are giving non-zero contri- 
butions and they can be presented as a sum of two terms, attributed to the different types 
of structural correlations in the system 

N 

r{q) = {zr Y.^l{q) + 2.Y,Ll,u,{q) , n = 21. (4) 

,k=l k2>k\ 



Here the first sum over ki = k2 = k contains the contributions L^{q) from N individual 
clusters, and the second sum over ki ^ ^2 contains the cross-terms L^^ ^.^ (g) corresponding 
to correlation of different clusters. In 2D systems these two terms L^{q) and fea^^) 
be presented in the following waj^ 



LUq) = j y"rfrirfr2Pfe(ri)pfe(r2)J„(g|r2i|)e-*"<^--, (5a) 




Ll,kM= dr,dr2PkA^i)pk2ir2)Jn{qmik,\)e "'=-^1. (5b) 



Here fc^ = Ra;2,A:i + ^21, where Rfea.fci is the 2D vector connecting centers of different 
clusters, ri and r2 are the 2D vectors defined in the local coordinate system (with its origin 



located at the center of k-th cluster). In Eqs. (5a, 5b) ^roi and 0^21 are the angles of 

1_ I I I k2,ki 

the vectors r2i = r2 — ri and Ullk^ in the external coordinate system, pfc(r) is a projected 
electron density of a cluster, Jn{qr) is the Bessel function of the first kind of integer order 
n, and the integration is performed over the area of a cluster. 

If all clusters in the system have the same internal structure but are randomly oriented 



in 2D space, the phase (f)r2i in the exponent of Eq. (5a) can be defined as 0r2i = 4>k + 0r2ii 
where (pk is the rotation angle of the fc-th cluster with respect to the fixed (reference) angular 
orientation 0°^^ of the clustei'^Sl. In this case, the term L^{q) can be expressed as follows 

Lm = e-"^'=L"(g), (6) 

and the integral L'^{q) is the same for all clusters: 

i^"(g) = J y"rfrirfr2p(ri)p(r2)J.(g|r2i|)e-^"<^°2i. (7) 

It is well seen from the structure of the integral in Eq. ([t]) that there is a direct cor- 
respondence between the electron density p(r) of LSs and the values of i^"(g). If L"-{q) is 
known, certain information about the electron density p(r), for example symmetry of clus- 
ters, can be obtained. Importantly, the first term in decomposition ^ does not depend 



on inter-particle distance and will have the same contribution for dilute and dense systems. 



Contrary to that, according to (5b) the term L^^^^(g) contains the inter-particle distance 
R'fc2,fci cLS an important parameter. Clearly, the contribution of this term will be different for 
dilute and dense systems. 

In the case of a dilute disordered system the contribution of the second sum in Eq. Q 
can be neglectecP^, and substituting Eqs. Q and ^ in Eq. ([s]) we obtain for the Fourier 
components of the CCF 

C-,=C-{q)-Al (8) 



where 



C'{q) = \L-{q)/L\q)\\ (9) 



is the structural term, which is directly related to the structure of a LS, and is the 
statistical term, which depends on the orientational distribution of LSs in the system. It is 
defined as an amplitude of the random phasor sumP'^, 

1 ^ 

A„ = A„e^^" = -^e^"'^^ (10) 

k=l 

where 9n is the phase of this sum. It is clear from Eq. (|8| that the value of the Fourier 
component in dilute systems is determined by the product of the structural and statistical 
terms. 

We analyzed the statistical behavior of for different angular distributions of cluster 
orientations in our previous papei'^. Our analysis has shown that in the case of a completely 
oriented dilute system (all (pk = <Po)i the statistical term is equal to unity, and the Fourier 
components reach their maximum values, determined by the structural term = C^{q)- 
If the distribution of orientations of LSs is Gaussian, with the standard deviation and a 
zero mean (0) = 0, the average value of the statistical term is 

{AD = exp (-nVJ) (1 - l/N) + l/N. (11) 



It follows from Eq. (11), that in the limit — t- oo, the Fourier components of the averaged 



angular CCF are defined by [see for details section VI] 

{C-) = C-{q)e^^{-n'al). (12) 
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According to these results for a dilute system with Gaussian distribution of orientations, 
the contribution of higher order Fourier components (corresponding to large n-values) will 
be much lower compared to the contribution of the lower order Fourier components. 

In the case of a uniform distribution of a large number of orientations, the statistical 
term fluctuates around its mean value (A^) with the standard deviation ctaz: 

{Al) = l/N, aAi = l/N. (13) 

It means that in this case fluctuations are of the same order as the mean value and their 



ratio does not depend on the number of clusters in the system. According to Eqs. (13) and 



(|8]), in the limit of a large number of orientations (A^ — )■ oo) the Fourier components of the 
averaged CCF have vanishing values, (C^ ) — )■ 0. 

In the case of a dense disordered system, when the average distance between LSs is of 
the order of the size of a single cluster, the second sum in Eq. ^ can not be neglected. 
Importantly, for such a dense system it can significantly affect the Fourier spectrum of the 
angular CCF. Taking into account both terms of Eq. (|4]), the Fourier components of the 
angular CCF (g can be written as the following sumP^l 



r\q) 



where 



SJ} + S^= 4 • Re 



k2>ki 



and 



(14) 

(15a) 
(15b) 

(15c) 
(16) 



NL\q) + 2Y,Ll^kM ■ 

k2>ki 

The summation in Eqs. (15a 15c) is performed over all LSs. In Eq. (14) the term 5" is 
defined by the local structure, and terms 5*2 , and contain contributions from the inter- 
particle correlations due to the second term in Eq. (|4]). Below we will explicitly investigate 
contribution of these terms to the total value of the Fourier components of the CCF. 



III. MODEL OF A DISORDERED SYSTEM AND SCATTERING GEOMETRY 



(a) Model 



The model and the geometry of the coherent x-ray scattering experiment used in our 
simulation^ is shown in Fig. [l|a). A coherent x-ray beam with a wavelength A = 1.5 A 
kinematically scatters from a 2D disordered sample of the total area 10 x 10 fim^. Diffraction 
patterns are recorded on a 2D detector of the size D = 12 mm (with the pixel size p = 
10 /im), positioned in the far-field at the distance L = 2.5 m from the sample. The chosen 
scattering geometry allows to cover a g-range up to 0.1 nm^^ (that corresponds the maximum 
scattering angle, 6'max = 0.14°), and also corresponds to about four pixels per speckle in a 
diffraction pattern for a given sample sizd^. 

As a model of a disordered sample a 2D system composed of identical clusters, generally in 
random positions and orientations, is considered [Fig. [TJ^b)]. We chose a centered pentagonal 
2D cluster with a 5-fold rotational symmetry for our simulations [see the inset in Fig. [l|b)]. 
This symmetry belongs to the forbidden motifs in the long-range crystalline matter but may 
exist on the short length scales in disordered system^, and is, clearly, most intriguing for 
a detailed study. Each cluster is composed of close-packed colloidal spheres with a radius of 
100 nm giving the total cluster size d = 600 nm. The sample is oriented with its 2D plane 
perpendicular to the direction of the incident x-ray beam [see Fig. [l]^a)], and the 5-fold 
rotational axis of each pentagonal cluster in the sample is parallel to the direction of the 
incident x-ray beam. The choice of parameters for our model was selected in analogy to the 
system studied experimentally in Ref.'^. 

In our simulations disordered samples are generated in the following way. Starting with 
an ordered 2D lattice of clusters, positional disorder is achieved by a sequence of random 
movements of all clusters. A single randomization step consists of Voronoi tessellatiorP^ 
of the sample area, with each Voronoi cell containing only one cluster, and a subsequent 
displacement of each cluster by a random vector within the Voronoi cell. After 10 — 50 
steps there is no positional order in the system. This procedure is repeated several times to 
generate the sequence of uncorrelated systems. Orientational disorder of clusters is imple- 
mented on the final step by applying a certain angular distribution of clusters to each of the 
generated samples. In Fig. [l]^b) one of the typical realizations of this procedure is presented. 
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(b) Structural term 

As discussed earlier, the contribution of a single cluster to different Fourier components of 
the CCF is determined by the term C'^{q) in Eq. (|8]) and strongly depends on its symmetry 
and structure. In case of a centered pentagonal cluster considered in our simulations due to 
its 5-fold symmetry only Fourier components with n = 101, where / is an integer number, 
have a non-zero contribution (see Appendix A). Importantly, the same term in Eq. ([s]) 
determines the g-dependence of Fourier components of the CCF. In Fig. [2] the values 
of |L"(g)p and |L°(g)p normalized by the form- factor of a sphere \f{q)\'^ as well as their 
ratio C^{q) = |L"'(g)/L°(g)p as a function of q are presented for n = 10,20 and 30. As a 
general rule, at low g— values, due to the properties of the Bessel functions, the strongest 
contribution to the CCF comes from the Fourier component with the lowest n {n = 10 in our 
case). At higher g— values contribution of other components becomes comparable. Notice, 
that at some g— values the function C^{q) contain sharp maxima corresponding to g— values 
where \L^{q)\'^ has small values. In a dilute system where we can neglect inter-particle 
correlations the factor NLP{q) coincides with the structure factor S{q) in SAXS. It means 
that one can expect an enhancement of the structural term for small S{q) values. Similar 
observations were made in experiment!^. Examination of Fig. |2]^c) also shows that small 
variations of g— values can result in large variations in the values of C"'{q) and consequently 
in Cg. This is also similar to results obtained in Ref.'^. It can be also seen in Fig. that 
there are some regions in g where the contributions of all components are close to zero. We 



may expect, that at these specific g— values the contribution of other terms in Eq. (14), that 
are responsible for inter-cluster correlations in dense systems, can become significant. 

The values of the structural term £"(g) at five selected values of g [shown by arrows in 
Fig. |2]^c)] are presented in Table [ij In the following we will compare these values with the 
values obtained for different systems by means of the XCCA. 



IV. LARGE NUMBER OF CLUSTERS IN THE SYSTEM 

In this section we consider model systems with different characteristics of orientational 
disorder. We are specifically interested here in analyzing the case of a dense system, contrary 
to the analysis of the dilute systems performed in our previous papei'^. In the following 
we consider three systems consisting of 121 pentagonal cluster^ with the same positions 
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of clusters in all three systems but with different orientational disorder in each system [see 
Fig. [l](b)]. These three systems are: a completely oriented system of clusters (Fig. |3]), a 
partially oriented system described by the Gaussian distribution of orientations with the 
standard deviation cr^ = 0.05 ■ 360°/5 = 3.6° and zero mean (Fig. [4]) and one with a 
uniform distribution of orientations in the angular rang^^fl] (— vr/S, 7r/5) (Fig. [s]). The relative 
separation between the clusters in all three cases is characterized by the ratio {R)/d ~ 1.5, 
where (R) is the average distance between the clusters estimated from {R) ~ 1/ ^/r^), and 
hq is the number of clusters per unit area (number density). 

Diffraction patterns for each case are presented in Figs. |3](a), |4]^a) andjsj^a). The most 
prominent features for all diffraction patterns are concentric rings appearing due to the 
contribution of the form-factor /(g) of a single colloidal sphere to the scattered intensity. The 
speckles superimposed on the diffraction patterns originate from the interference of x-rays 
coherently scattered on the spatially disordered clusters, distributed within a finite sample 
area. The corresponding angular averaged intensities {I{q, V^))<^ for each diffraction pattern, 
as well as the angular distribution of orientations of clusters in each system are presented 
in Figs, [sj^b), [4]^b) andjsj^b). The following momentum transfer g— values: qi = 0.023 nm~^, 
q2 = 0.029 nm~^, q^ = 0.036 nm~^, q^ = 0.043 nm~^ and gs = 0.059 nm~^ were considered 
for the detailed analysis in each case [shown in Figs. [3]^b), |4](b) andjsj^b) by red circles]. The 
CCFs Cq(A) as well as corresponding Fourier components were evaluated using Eqs. ([T| 
2b) for each system at these selected g— values and are presented in Figs. |3](c,d), |4](c,d) and 
gc,d). 

Comparison of all three cases shows that similar to a dilute system the dominant Fourier 
components {n = 10 and 20 in our case) are well defined for more oriented systems (Fig. [3] 
and Fig. |4]). It is also well seen that for a Gaussian distribution of orientations (Fig. |4]) 
at higher g— values (g = gs in this case) the magnitudes of the Fourier components with 
n = 10 and n = 20 are significantly reduced in comparison to a completely oriented system 
(Fig. [3]). This decay can be attributed to the exponential factor present in the statistical 
term in Eqs. ( Tlp2 ). 



Our analysis of the first two systems with the high degree of orientational order shows 
that the symmetry of clusters composing the system can be well determined by means of 
the XCCA analysis. However, the absolute values of Fourier components can be strongly 
affected by the inter-particle correlations (see Tables [H[ll|). Comparison of the absolute values 
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of the Fourier components at n = 10, 20 for the completely oriented system considered 
here shows that they may fluctuate around the value of the structural term C^{q) for about 
few tens of percent at different g— values (see Table [l]). Results of similar calculations for a 
dense system with the Gaussian distribution of orientations are presented in Table [TT} They 
are compared with the averaged Fourier components of CCFs defined in Eq. ( [l2| . They 
show similar fluctuations of the values of due to inter-particle correlations in the system 
and statistical fluctuations of orientational distribution of clusters (see next section for a 
detailed discussion). 

Our simulations for the systems with different positional disorder have shown that 
also depend on the particular realization of the system. Such behaviour of can be 
explained by the specific constructive and destructive interference of the scattered x-rays in 



each particular system. Due to the interference term (5*2 + S'^) in Eq. (14) the inter-particle 
contribution can increase or decrease the contribution of the structural term S*" in a dense 
system. Note, that the term 5*4 is always non-negative and can only increase the value 
of Cq. That brings us to the conclusion that determination of the absolute values of the 
Fourier coefficients C^, or |/"'(5')| that is necessary for the direct structural determination 
of the clusters by phase retrievaP can be strongly affected in close packed systems. 

Contrary to these two cases, the system with the uniform distribution of orientations 
shows quite different behavior. As a result of the complete orientational disorder, the CCFs 
shown in Fig. [5|^c) do not represent a simple ~ cos(nA) form but have rather complicated 
angular dependence that can be also well observed in the corresponding Fourier spectra in 
Fig. |5]^d). For example, the component with n = 8 (and some other higher values of n) 
dominates at g = gi contrary to the previous more ordered systems, where the component 
with 72 = 10 was strongly dominant at that g— value. It is well seen in Fig. |5]^d) that in 
this case there are many Fourier components with comparable values in a wide range of the 
spectrum. We should note also here that the absolute values of all Fourier components in 
this case are about one order of magnitude lower then in the case of the ordered systems [see 
Figs. [3]^d) and|4]^d)]. This general behavior of disordered systems can be also explained by 
examining the contribution of two terms in Eq. (|4]). In the system with the uniform distri- 
bution of angular orientations the absolute value of the first term in Eq. (|4]) is significantly 



reduced due to the asymptotic behavior of the statistical term given in Eq. (13). For 



a large number of orientations it drops as 1/A^ and becomes much smaller then the second 
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term in Eq. Q that is determined mostly by the inter-particle correlations. In this case 
inter-particle correlations give the major contribution to CCFs and determination of the 
structure or symmetry of the clusters composing the system becomes practically impossible. 

The results of our simulations for the 2D samples with different degree of orientational 
disorder presented in this section clearly indicate, that the possibility to extract the infor- 
mation on the symmetry of individual clusters in the systems with a large number of clusters 
by the cross-correlation analysis is strongly related to the degree of orientational disorder in 
the system. The preferred situation is when the system is partially aligned, while correlation 
analysis is practically impossible in a completely disordered dense system. 



V. ENSEMBLE-DEPENDENT FLUCTUATIONS OF CORRELATIONS 



In this section we will consider in more detail the statistical behavior of the CCFs as a 
function of orientational disorder in a system with a small number of possible orientations. 
As was discussed above, in a dilute disordered system composed of a certain number of 
clusters the statistical term in the CCF is fluctuating around an average value (Af^) 



[see Eqs. (11) and (13)], determined by the statistics (Gaussian, uniform, etc.) and number 
of angular orientations of clusters. Contrary to the previous section where we analyzed 
systems with a large number of clusters (A^ = 121), here we consider systems with much 
smaller number of clusters for which statistical fluctuations of the scattered intensity become 
considerably larger. 

We analyze here three samples, each composed of = 11 clusters. The positions of 
clusters are the same in each sample, and their relative separation is given by the ratio 
{R)/d ~ 5.0. We consider for each sample a uniform distribution of 11 orientations of 
clusters within the angular range (— vr/S, tt/S) [see Fig. [6]. 

The results of the Fourier analysis of the CCFs for these three systems at the selected 
values of q are shown in Fig. |6] (note that a different vertical scale for the Fourier spectra was 
used at different q values). Comparing the Fourier spectra calculated at the same q value for 
these systems one can see, that they significantly vary only due to a different orientational 
distribution of clusters from one sample to another. In the case of dilute systems, fluctuations 
of the values of Fourier components are defined by the statistical distribution of orientations 
of clusters [see Eqs. ([8 13)]. In the systems considered here, these fluctuations are defined 
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by the interplay of two factors, the distribution of cluster orientations [term Si in expansion 



(14)], and their relative spatial positions (terms S'2, 5*3 and 5*4). This leads in some cases 



to a large contribution of S"" , sufficient to dominate over the other terms in Eq. (14) [see, 
for example, the results for the Sample 3 in Fig. [6]^c)], and in some cases to a smaller 
contribution, allowing the inter-cluster correlations to dominate in the Fourier component 
Cg [see, for example, the results for the Sample 1 in Fig. |6](a)]. Therefore, even if the number 
of possible orientations of clusters in a system is small, it may be not possible to extract 
the information on local structure of such a system from a single diffraction pattern. The 
statistical behavior of CCFs discussed here could give a possible explanation to experimental 
observations in Ref.'^, where diffraction patterns taken at different times show different types 
of angular correlations. 

As demonstrated here in a system with a small number of randomly oriented clusters 
reliable information on the local structure and symmetry of the clusters can not be obtained 
by the analysis of a single diffraction pattern due to statistical fluctuations of intensity. We 
discuss in the next section how this problem can be solved in certain cases by averaging of 
CCFs. 



VI. EFFECTS OF ORIENTATIONAL DISORDER ON THE AVERAGED CCF 

In our previous discussion we presented the results of the XCCA applied to the analysis 
of the diffraction data obtained for a single realization of a disordered system. We showed 
that statistical fluctuations of CCFs from one realization to another is a strong limitation for 
the determination of reliable information on cluster symmetry and structure. To overcome 
this difficulty it was proposed, first by KamPS, to analyze CCFs averaged over an ensemble 
of diffraction patterns, instead of analyzing CCFs calculated for a single diffraction pattern 
(see also recent publication^^'^'^ . Below we will investigate conditions when this approach 
is valid. 

We will consider here the CCF (Cg(A))j\,/ averaged over a sufficiently large number M of 
diffraction patterns, 

M 

n(A))M = l/M^C7(A), (17) 

m=l 

where C^(A) is the CCF obtained for the m-th diffraction pattern [Eq. ([T|]. Performing the 
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Fourier transform of both parts of Eq. (17) and taking into account that Fourier transform 



is a hnear operator we obtain for the Fourier components of the averaged CCFs 

M 

n")M = i/Mj]{Qr, 

m=l 



(18) 



where {C"}™ are the n—th. Fourier components of the m-th CCF C™(A). Eq. (18) means, 



that the n—th Fourier component of the averaged CCF is an average of the corresponding 
Fourier components determined for each diffraction pattern. Taking into account expression 



(|14|) for each individual Fourier component of the CCF we obtain 



K) 



M 



Piq) 



M 



+ 



M 



\P{qW 



M 



(|/o(g)P 



on 
O4 



(19) 



M 



where S"' are defined by Eqs. (15a 15c) and |/°(g)p by Eq. (16). Below we show by direct 



simulations that for a sufficiently large number of diffraction patterns M the ensemble 
averaged term ((5*2 + •S'g )/|/°(g)p)M asymptotically approaches zero and we have for the 
Fourier components of the averaged CCF 



where 



(dust) (int—clust) 

(dust) Ml 
{Cg (int-clust)) M = {S2/\I^{q)\'^)M- 



(20) 

(21a) 
(21b) 



This result indicates that in a general case for a sufficiently large number of diffraction 
patterns the Fourier components of the averaged CCF can be represented as an additive sum 
of two positive valued contributions. The first one is given by the structure of individual 
clusters and the second one is determined by the inter-cluster contribution. 

In the case of a dilute disordered system, when the inter-cluster contribution {C"^ (^^^t-ciust)) m 



in Eq. (20) can be neglected, we obtain 



{Cq)M - {Cq{clust))M - '^"(g) ■ (A^)m, 



(22) 



where {A^)m is an ensemble averaged square amplitude of the random phasor sum (10) 



M 



{aI)m = i/mY,{aIY 



(23) 



m=l 
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For a sufficiently large number M of diffraction patterns, the value of (A^) m approaches 
its statistical limit (A^) [see Eqs. (11) and (13)]. This result immediately explains the ad- 



vantage of calculation of the average CCFs in dilute systems as soon as it converges to a 
finite not fluctuating value. This opens the opportunity to determine the structural term 
for systems with orientational disorder. At the same time our analysis shows that 
the favorable situation for such averaged analysis are dilute systems with a small number of 
orientations. If the number of orientations increases, contribution of this averaged statis- 
tical term {ADm decreases [for example, as ~ 1/A^ for a uniform distribution of orientations 



according to Eq. (13)], and the structural contribution {C'^ (^lust)) m can become smaller than 
the inter-cluster contribution {C^ [i^t-ciust)) m ■ In Appendix |b] we give a detailed derivation 
for the statistical estimate of the inter-cluster contribution to {Cg)M [Eq. (20) . 

One important question that is relevant to the discussion in this chapter is how many 
diffraction patterns M will be sufficient for calculation of the average (A^) m in order to 
reach reliably its statistical limit {A"^). In other words, what is an averaging accuracy e 
for a number M of diffiaction patterns. This question can be addressed by establishing a 
confidence interval for {A'^)m- Our analysis has shown (see Appendix [C] for details) that 
in the case of a uniform distribution of cluster orientations and averaging accuracy of 10% 
the probability that {A'^)m lies in the interval from 0.9(A^) to l.l(A^) is more than 90% for 
M = 10^, and more than 99% for M = 10^ diffiaction patterns. 

We demonstrate these results by analyzing the ensemble averaged CCFs for three systems 
consisting of (a) 11, (b) 60 and (c) 121 clusters with the relative distances (a) {R)/d ~ 5.0, 
(b) {R)/d ~ 2.2 and (c) {R)/d ^ 1.5 respectively. For each realization of a system in 
the ensemble different set of angles and cluster positions were considered. In Fig. [7] the 
results of the Fourier analysis of the CCFs averaged over M = 1000 diffiaction patterns 
are presented. According to our previous analysis, for this number of diffraction patterns 



the term (Cg (ci„<ij))M in Eq. (20) with the probability more then 90% lies in the interval 

Comparison of the averaged Fourier spectra for the system containing 11 clusters 
[Fig. |7|a)], with the ones calculated for different single realizations of the same system 
[Fig. [6] shows a strong enhancement of the contrast of the contribution from the internal 
structure of the pentagonal cluster, as compared to the inter-cluster contribution. At the 
same time, comparison of the magnitudes of the same Fourier components with the ones 
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obtained in a dilute limit for a system with the uniform distribution of cluster orientations 



[see Eq. (22) with M — )■ cxd] reveals significant deviations [see Table III . We attribute this 
effect to the presence of the inter-cluster contribution in the latter case. 

Comparison of the results presented in Fig. [7] for three different samples shows, that 
as soon as the number of clusters in a dense system increases, the Fourier components 
{dust)) M related to the local structure decrease following the dependence [see 

Eq. (13)]. It is well seen in Figs. [7|^b) and [7](c) that for a system with the number of 
clusters = 60 the Fourier components with n = 10 and 20 are already hardly resolved 
over the inter-cluster contribution, and for the system with = 121 clusters their contribu- 
tion is not resolved over the level of inter-cluster contribution. Therefore, if the number of 
possible orientations of LSs in a dense system is sufficiently large, as in Fig. [Tj^c), calculation 
of the average CCFs may not provide information on the local structure in a disordered 
system. 

The contribution of the inter-cluster correlations {Cq i^int-dust)) averaged CCF 
is well seen in an entire Fourier spectrum for all three systems presented in Fig. [7} In 



Appendix [B] an asymptotic estimate of the inter-cluster contribution in {Cq)M [Eq- (20)] is 
obtained (see solid red line in Fig. [?]). As one can see in Fig. [7| our estimate of the inter- 
cluster contribution quite accurately reproduces the results of the direct calculations for all 
systems considered here. Our analysis has shown that for a sufficiently large number M of 
diffraction patterns the magnitude of the averaged inter-cluster contribution {C'^t^^^t-dust)) ^ 
depends on the shape and size of the sample, as well as on the size of particles, and rather 
weakly on their number A^. 

The evolution of different terms in the expansion of the averaged Fourier component 



{Cq)M [Eq. (19)] as a function of the number M of diffraction patterns is presented in 
Fig. |9] The results are shown for with n = 10 for the same systems and q values as 
in Fig. [7| These results demonstrate the convergence of different terms to their average 
values. In particular, as one can see from Fig. ^ the averaged term ((S'2 + S^)/\P{q)\^) m in 



Eq. ( 19 ) asymptotically reaches zero and can be neglected after averaging over a few hundred 
diffraction patterns. For a system with the number of particles A^ = 11 the structural term 
strongly dominates over the inter-cluster contribution at most of the q values (it is only 
lower for q = qs). Contrary to that case, for the systems with the number of clusters A^ = 60 
and A^ = 121 the inter-cluster contribution term is larger than the term corresponding to 
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the structure of individual clusters for most of the g— values. That is the reason why the 
structural contribution could be hardly resolved for a system with N = 121 clusters in 
Fig. Be). 

Our results clearly show that, in the general case, the ensemble averaged CCF {Cij{A))m 
contains both the contribution from the individual clusters and the inter-cluster contribution 



[see Eq. (20)]. Therefore, the ensemble averaged CCF {Cq{A))M determined from a system of 
coherently illuminated clusters, in general, can not be considered as a function of only 'single- 
particle quantities' as it was stated in Refs.'^'^. For a large number N of cluster orientations 
in the system the inter-cluster contribution {Cqf^int-dust))^ '^^^ become dominant over the 
structural contribution {Cg(^ciust))^i- next section we consider a partially coherent 

scattering from disordered systems as one of the possible ways to suppress the inter-particle 
contribution to the CCF. 



VII. EFFECTS OF PARTIAL COHERENCE 

In the previous sections we analyzed diffraction data obtained under conditions of coher- 
ent illumination of a disordered 2D sample. The influence of partial coherence of x-rays on 
the results of the x-ray cross-correlation analysis was briefly discussed in Ref.'^. Here we 
will consider this question in more detail. 

Partial coherence of x-rays can be characterized in terms of the mutual intensity 
functionl^^'Sl J-^^(^ri,r2) of the beam incoming on the sample and defined as Jin(ri,r2) = 
{E{ri,t)E{r2,t))T, where averaging is performed over times T much longer than the fluctu- 
ation time of the incoming x-ray field. It describes the statistical properties of the wavefield 
as a correlation function between two values of the electric field, E{ri,t) and E{r2,t), at 
different points ri and r2 in space and at the same time t, averaged over fluctuations of the 
wavefield. It is convenient also to introduce a normalized complex coherence factor, 

/iin(ri, rs) = Jin(ri, r2)/[/in(ri)/in(r2)]'/', (24) 

where /in(i'i) and /in(i"2) are the intensity values of the incoming beam at points ri and r2. 



respectively, averaged over fluctuations of the wavefield. Using the definition (24), the effect 
of partial coherence of x-rays on the distribution of the scattered intensity /pcoh(q) can be 
expressed in the far-field limilP^as a convolution of the coherently scattered intensity /coh(q) 
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with the Fourier transform /iin(q) of the complex coherence factoi'^, 

Voh(q) = /coh(q) ® Aiin(q). (25) 

One can directly use this equation for calculations of the distribution of the scattered inten- 
sity for partially coherent illumination of a sample, as soon as coherence properties of the 
incoming beam are defined by the known complex coherence factor. 

We will demonstrate the influence of partial coherence of x-rays on the results of the 
XCCA by analyzing the scattering experiment with different degree of spatial coherence of 
the incident x-rays for a disordered system consisting of 11 clusters, presented in Fig. |6](c) 
(for Sample 3). We assume in our simulations a Gaussian form of the complex coherence 
factor (so-called Gaussian Shell-modeP^ 

/iin(ri,r2) = exp[-(ri - v^f/'^llo^.l (26) 

where Zcoh is the transverse coherence length. Three different values of the transverse coher- 
ence length, /coh = 1-2 /im, 600 nm, and 300 nm are considered in our calculations. These 
values of the coherence length /coh allow to probe typical length scales in a chosen sample, 
starting from comparably large coherent length, going down to the size of a single pentago- 
nal cluster and below. In Fig. |8] diffraction patterns calculated according to Eq. 
a function of the transverse coherence length /coh, and corresponding Fourier spectra of the 
CCFs calculated at three different q values for each pattern are presented. As one can see 
from Fig. |8j decrease of the transverse coherence length corresponds to smearing of speckles 
in a diffraction pattern. At the same time, smearing of speckles reduces the contribution 
from inter-cluster correlations in the Fourier spectra of the CCFs. Partial coherence of the 
incident x-rays acts as a filter for the Fourier components of the CCF with a threshold de- 
fined by the value of the transverse coherence length. Decrease of the value of /coh down to 
the size of a single pentagonal cluster considerably reduces the high-frequency contribution 
in the Fourier spectra, however, the contributions from the inter-cluster correlations do not 
vanish completely [compare Fig. |8](b) with the results for the case of fully coherent illumi- 
nation of the Sample 3 in Fig. [6]^c)]. Further decrease of the value of /coh below the size of a 
single cluster leads to a further refinement of the Fourier spectra, and the contribution from 
the local structure becomes clearly dominant [Fig. [8|c)]. 

Comparison of the Fourier components of the CCFs {n = 10, 20) for the case of 
fully coherent scattering (C^ /coh = oo) and partially coherent scattering (C^ /coh = 
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25p6) 



as 



600 nm) from a system of 11 clusters with a uniform distribution of orientations is presented 
Note, that absolute values of the Fourier components corresponding to 



in Table 



IV 



partially coherent data are significantly lower at all q values, as compared to '^"^ values. 
At the same time, the inter-cluster contribution is suppressed significantly stronger and as 
a result one gets higher contrast for the LS contribution. 

It is interesting to compare these results with the ones presented in Ref.l^S, where the 
effects of a partial coherence were studied on a model sample represented by an ordered 
array of identical particles in the same orientation. It was demonstrated^^, that when the 
transverse coherence length of the incoming beam illuminating such a sample approaches the 
size of a single particle, the resulting diffraction pattern looks similar to the one produced 
by a single particle. Taking this into account, one would expect that for the disordered 
sample, considered here, a decrease of the transverse coherence length down to the size 
of a single pentagonal cluster will lead to a diffraction pattern, which resembles a sum of 
11 diffraction patterns obtained for a single cluster in 11 different orientations. Such a 
diffraction pattern with the partial illumination of x-rays could be sufficient for the future 
phase retrieval analysi^ ^ * ^ * ^ * ^ l 

Analysis of the effect of partial coherence on the results of the cross-correlation analysis 
allows the following conclusions to be made for disordered 2D systems. If the value of the 
transverse coherence length is of the order of a size of a single particle (cluster, molecule), 
the contribution of the inter-particle correlations strongly decreases (completely vanishing in 
the case of a dilute system), enabling easier identification of the Fourier components related 
to the internal structure of particles. Therefore, one may prefer to perform a scattering 
experiment with partially coherent x-ray beam, in order to extract the information on local 
structure of a disordered system. On the other hand, varying a degree of coherence of x-rays, 
one could probe the structural correlations in a disordered system at different length scales, 
which may be particularly useful in studies of a medium range order. 



VIII. CONCLUSIONS 

In a summary, our simulations have demonstrated that in order to get reliable information 
about the local structures the systems with a high degree of orientational order are preferable 
ones comparing to the completely disordered systems. Due to the statistical nature of the 
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CCFs the information on local structure of a completely disordered system from a single 
diffraction pattern is practically not accessible. Averaging of CCFs over an ensemble of 
diffraction patterns allows to overcome this difficulty and gives reliable information on the 
structure of clusters forming the system. An estimate of a number of diffraction patterns 
necessary to determine an averaged CCF with a certain accuracy was obtained. Our analysis 
shows that few tens of thousands of diffraction patterns have to be averaged to get a reliable 
result. 

The limiting factor for this analysis becomes the number of particles in the system and 
the inter-particle correlations that are always present in the conditions of coherent illumina- 
tion. Our simulations show that the systems with a small number of clusters are preferable. 
Averaging of CCFs over an ensemble of diffraction patterns for such systems can signifi- 
cantly enhance the contrast of the structural contribution as compared to the inter-particle 
contribution. For the systems with a large number of clusters this averaging procedure 
could become not efficient due to the strongly suppressed signal (it scales as with the 
number of particles A^) from the structural contribution. Our analysis has shown that the 
averaged inter-particle contribution practically does not depend on the number of particles 
and their density, but rather depends on the size and the shape of the system as well as on 
the size of clusters composing the system. Utilizing partially coherent incident beams with 
the coherence length about the size of the clusters composing the system could suppress 
the inter-particle contribution and make possible determination of the local structure of a 
disordered system. 

Interesting open question left for the future work is an extension of XCCA to the three- 
dimensional (3D) systems. The additional degree of freedom for orientation of 3D clusters 
introduces a significant complication in the extraction of the structural information from 
the XCCA analysis. As it was pointed out in our previous workp! one possible solution can 
be the measurements performed at high scattering angles when the curvature of the Ewald 
sphere has to be taken into account. This will be especially interesting case for atomic 
systems (organic and non-organic molecules, biological systems, etc.) as opposed to the 
systems composed of colloidal particles considered in this work. The full cross-correlation 
analysis for such practically important systems could open new horizons in exploring the 
nature of disordered systems. 
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Appendix A: 



Here we present the results of calculations of the function L"'{q) [Eq. ([7])] for a centered 
pentagonal 2D cluster [Fig. [l|^b)]. The direction of the incident x-ray beam is assumed to 
be parallel to 5-fold rotational axis of a cluster. We define the electron density of a cluster 
as a real-valued quantity in the following form 



(Al) 



1=1 



where fi{q) and are the form factor and the radius vector of the i-th scatterer and Ng is 
the number of scatterers in the cluster. Substituting Eq. (Al) into Eq. ([T]) and assuming 
that the cluster is composed of scatterers with the same form factor f{q) we obtain 

L-{q) = \m\'J2Uq\rst\)e-'-'''^'^ (A2) 



s.t 



Using the coordinates of scatterers in Eq. (A2) we obtain for a centered pentagonal 2D 
cluster 

' |/(g)P{6J„(0) + 10[J„(ga) + J„(Aga) + J„(A2ga)]} if n = 0, 

n mod 20 = 0; 

= <! 10|/(g)|2[J„ (qa) - J„ {A^qa) - J„ (Asga)] if n mod 10 = 0, (A3) 

n mod 20 ^ 0; 

other n, 

where a is a distance between centers of the nearest particles inside the pentagonal cluster, 
Al = ^Jl{5 - V5) ^ 1.18, and A2 = y^^(5 + ^5) ^ 1.9. Note, that the function 
presented here for a centered pentagonal cluster differs from the one defined in Ref.I^ for a 
LS with a fivefold symmetry due to the missing central particle in the latter case. 
In our simulations we consider the following form factor of a spherical particle 

47r[sin(gr) — (gr) cos(gr)] 



fiq) = P- 



(A4) 



where p is the average electron density inside a particle and r is the radius of a particle. 



Appendix B: 



Here we determine a statistical estimate of the average Fourier components of the inter- 
cluster contribution {C^ (^int-dust)) i^^^ -^l- (20)]. First, we will consider a simphfied model of 
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a disordered system, which consist of spherical particles. Next, we will show by simulations 
that asymptotic equations derived for this simplified model can be used to determine the 
inter-particle contribution in a disordered system composed of more complicated objects. 

can be written as 



According to Eq. (21b) the normalized Fourier component i^-^t-dust) 
follows 



a 



q [int—clust) 



S2l[i\<i)f 



(Bl) 



where 



S2 = A 



k2>ki 



k\,k2 



NL\q) + 2 Ll^kM) 



k2>k\ 



(B2a) 
(B2b) 



Here was decomposed into contribution of two terms with ki = k2 and ki ^ k2 similar 



to the expression in Eq. Q, L^_^ f^^{q) is defined by Eq. (5b), and by Eq. with 

n = 0. 



Using the expression for the electron density p(r) [see Eq. (Al)] in Eqs. (5b) and ([T]) for 
a system of clusters composed of identical scatterers we can writd^ 



s,t 

L\q) = J2Mq\rst\). 



(B3a) 
(B3b) 



s.t 



In Eq. (B3a) the summation is performed over positions of the scatterers inside the clusters 



ki and /c2, and in Eq. (B3b) inside a single clusteilSS. 



Now, we consider a simplified model of a disordered system, where we substitute a cluster 
by a spherical particle of the same size as the original cluster. In this case r^t = and 



Eqs. (B3a, B3b) reduce to 



= 1. 



fc2,fei 



(B4a) 
(B4b) 



Further, the double summation in 5*4 and [Eqs. (B2a) and (B2b)] can be substituted 

by a single sum over the same set of vectors. Preserving the number of contributions in the 
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original sum we can write 



5r = 4|5^J„(g|R,|)e-*"'^-^ 



fc=i 



N + 2j2Mq\^k\) 



k=i 



(B5a) 
(B5b) 



where Ni = N{N -l)/2. 

The sum in Eq. (B5a) is, in fact, a random phasor sumP^ESl (defined as 

1 



(B6) 



k=l 



where and On are the length and phase of the random phasor sum B„, and the amplitude 
and phase of each term are given by 



(B7) 



With these definitions we can also formally write the random sum in Eq. (B5b) in the 
following form, 

1 



k=l 



Substituting Eqs. (B6 - B8) into Eqs. (B5a, B5b) we obtain 



S2 = mfBi 

I 2 



(B9a) 
(B9b) 



=[N + 2N^BoY = + ANN^Bq + AN^B^ 
As soon as we are interested in averaged characteristics of the Fourier component 



(int-ciust) i^'i- (^1)] wc uotc that both terms 5*4 and [/°(g)]^ can be considered as 



random variables as well. Applying the theorem on a linearized approximatiorP^ to evaluate 
{^q{int-dust)) cxpress an average of the ratio of these two variables as a ratio of the 

corresponding averages 



^ {s2)/{{i\q)r), (BIO) 

where 

{S2)=4N',{Bl), (Blla) 

{[I\q)?) = N^ + iNN^iBo) + AN^{B^). (Bllb) 
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It is clear that this approach is vahd for g— values when ([J°(g)]^) 7^ 0. Our direct simulations 



have shown that approximation described by Eq. (BIO) is well satisfied for our systems. 
According to Eqs. (BIO), (Blla) and (Bllb), in order to estimate the value {C"^ (^^nt-dust)) 



we need to determine the statistical averages {B^) and {Bq). Following the approach of 
RefP, one can obtain the marginal probability density function for the random variable B^ 
(n 7^ 0), that reduces for a large number of particles to a normal distribution 



2o-2 

^ n 



exp 



B^ 



20-2 



Bl > 0, 



where 



Ar2 



k=l 



m?) _ {JMm 

2 2Ni 



(B12) 



(B13) 



Using the probability density function P{B^) the mean value of the random variable B^ and 
its variance can be determined for n 7^ as 



{Bl)= 2al = {4{q\m/Nu 
^h='^cTt,= [{Jl{q\m/Ni] 



(B14a) 
(B14b) 



According to Eqs. (B14a) and (B14b) the value B^ fiuctuates around the mean value (B^) 



with the standard deviation cr^a. According to Eq. (Blla), this also determines the statis- 



tistical behaviour of the term S^- Note, that similarly to the random phasor sum considered 
for LS contribution, averaging over sufficiently large number M of diffraction patterns leads 
to a decrease of the variance and the result of such averaging {B'^)m asymptotically 
approaches its statistical limit, i.e. {B^)m {B'^) (see Appendix [C|) . 
The average values (Bq) and {B^) can be directly determined as 



(Bo)- 



1 

^ k=i 



(B15a) 



iB'o)=j,j:Kbi) 



1 



.k=l 



[{j'o{q\m + iNi-l){Mq\m'] 



(B15b) 



In the derivation of Eq. (B15b) we assumed that for ki 7^ k2, (Jo(g|RfcJ) Jo(Q'|Rfc2l)) 
(Jo(g|RfcJ))(Jo(g|Rfe|))- 
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Using these results in Eqs. (Blla) and (Bllb) we obtain for the statistical averages {S2) 
and {[I^qW), 

{S2)=m{j'^{q\R\)), (B16a) 
mQ)?) = N' + 4N, [N{Mq\m + (Jo'(g|R|)) + {N, - l)(Jo(g|R|))^] • (B16b) 
For a large number of particles in the system ^ 1, A^^i ^ and substituting now 



Eqs. (B16a, B16b) into Eq. (BIO) we obtain for the statistical average 



(-*-^'-*)^ 1 + 2iV( Jo(g|R|)) + 2( j2(g|R|)) + m{Mq\R\W ' 
Our direct simulations for the systems with a large number of clusters considered here show, 
that the terms Jo(g|R|)), (Jo(g|R|)) as well as A^^(Jo(g|R|))^ in the denominator of 



(B17) 



Eq. (B17) are much smaller than unity and Eq. (B17) reduces to the following asymptotic 



expression 

(Q(»t-d«st))^2(j2(g|R|)). (B18) 
As we can see from this expression, for the systems considered here the inter-cluster contri- 



bution practically does not depend on the number of clusters N. Substituting Eqs. (BIS) 



and (22) in Eq. (20) for a sufficiently large number of diffraction pattern we obtain 



{C:) = r{q)-{Al)+2{4{q\m- 



(B19) 



Taking into account that for a uniform distribution of orientations {A'^) = 1/N [Eq. (13) 
we obtain the upper limit for the number of particles A^, for which the structural contri 



bution dominates over the inter-cluster contribution. The first term in Eq. (B19) strongly 
dominates, if the number of clusters A^ satisfies the following condition 



Ar< 



(B20) 



For the systems considered here this estimate shows that the number of particles in the 
system should be considerably lower than about a few hundreds, that is in a good agreement 
with our observations. 

The averages ( Jo(q'|R|)) and (J^(g|R|)) (for all n values) required to calculate {C"^ (int-dust)) 
can be determined using the radial distribution function (RDF) P{R) 



Rn 



P{R)MqR)dR, (J^(g|R|)) 



R„ 



P{R)Jl{qR)dR, (B21) 



R„ 
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where -Rmin and -Rmax are the minimum and maximum inter-particle distances in the system. 
According to its definition!^, P{R) gives a number of particles in an annulus of thickness 
dR at a distance R from another particle. Note, that P{R) defined in such a way should 



be normalized by the number of particles before using it in Eq. (B21), in order to satisfy 
the condition J^™'"' P{R)dR = 1. In our simulations we consider a statistical distribution 
of spherical particles with a diameter d within the square 2D sample with the length / of 
a side of the square. For such a system, it is convenient to determine the RDF using the 
"partial radial distribution function"H2l (PRDF) Ppart(-R, m, w), 

l>l-d/2 i>l~d/2 

P{R) = — / Ppart(i?, u, iy)dudu, d<R< V2{1 - d), (B22) 

^fs Ju=d/2 Ju=d/2 

where Ppart{R,u,i>) defines the number of particles in an annulus of thickness dR at a 
distance R from a particle at a position P which is located from the sides of the square 
at distances u and v [see Fig. 10 , and Sfs = {I — d)"^ — Tcd"^ is a part of the sample area 



accessible for Ppart(-R, m, i^)- Depending on the radius R and the position P, Ppart(-R, i^) 



can be proportional to the length of the circumference [Fig. 10 a)], the length of a single arc 
[Fig. [Io|(b)], or the total length of the few arcs [Fig. [Io|(c)]. Clearly, the RDF depends on 



the size and the shape of the particles and the sample. Eq. (B22) ensures that the minimum 



distance between the particles is d, and minimum distance from a particle center to the 
sample edge is d/2; it also defines the region of integration in Eq. (B21) with Rmin = d, 
Rmax = - d). The RDF P(P) calculated using Eq. (^B22| for a system of spherical 



particles with d = 600 nm and the size of the sample / = 10 /im is presented in Fig. [10j(d). 
As we can see from this figure, the maximum of this distribution lies between 4 /im and 
6 fim that corresponds approximately to the half size of the system. 

We compare now the results of an asymptotic estimate {C"^ (^int-dust)) inter-cluster 



contribution for a system composed of spherical particles using Eqs. (BIO), (B16a) and 



(B16b) with the results of direct calculations for different systems by averaging (^^^t-ciust) 
in Eq. (Bl) over M = 1000 realizations. In the latter case, the coordinates of particles in 



each realization were used in calculations of Eqs. (B2a) and (B2b). Three 2D disordered 



systems composed of (a) spherical particles [Fig. 10 e)], (b) triangular clusters [Fig. 10 'f)] 
and (c) centered pentagonal clusters [Fig. [Io|^g)] were considered in these direct calculations. 
Fourier components {C^ ^^nt-dust)) m and calculated using Eqs. (IbT]) and 



(BIO) for n = 4 and different q values for systems with A^ = 11, 60 and 121 particles are 
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presented in Fig. 11 The Fourier component with n = 4 was considered because for the 



chosen symmetry of clusters only the term 5*4 in Eq. ( 14 ) contributes to [terms S*" and 
(S'2 + 5*3) are equal to zero]. In all three cases direct calculations show that the inter- 
cluster contribution practically does not depend on the internal structure of particles. As 



one can see from Fig. [TT| the inter-cluster contribution can be accurately estimated using 
the asymptotic expression. The inter-cluster contribution calculated using the asymptotic 



expression [Eqs. (BIO), (B16a) and (B16b)] for three systems considered here for an entire 
Fourier spectrum is presented in Fig. [7] with the solid red line. As one can see, our estimate 
quite accurately reproduces the results of direct calculations for all systems considered here 
in a wide range of n values. We want to note here, that asymptotic values for the system 
with = 11 clusters slightly overestimate contribution from the inter-cluster correlation 
[see Fig. [7|^a) and Fig. 11 a)] that is due to the violation of the central limit theorem, used 



in the derivation of Eq. (BIO). 

Finally, our analysis shows that the asymptotic expression for the average inter-particle 



contribution {Cq {j^nt-dust)) [Eqs. (BIO), (B16a) and (B16b)] can be effectively used as an 
estimate of this contribution to the total Fourier component {Cg)M in Eq- (20). We demon- 
strated that {Cq (^int-dust)) depends on the shape and size of the sample, as well as the size 
of particles, and rather weakly on their number A^ in the system. 



Appendix C: 

Here, we give an estimate of the number of diffraction patterns M that is required to 



determine the ensemble averaged amplitude {ADm [Eq. (23)] with a given accuracy e. We 
define e as a maximum deviation of {A'^)m from its statistical limit {A"^). 

Each single x-ray pulse produces a diffraction pattern from a static distribution of clusters 
in the system, and for the m-th diffraction pattern the square modulus Af^{m) of the random 



phasor sum (10) can be written as 

Alim) = {Al)+SAl{m), (CI) 

where {A"^) is the statistical average and 5A^(m) describes fluctuations of this average value 
from one diffraction pattern to another. 
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Introducing an average over M diffraction patterns one can write 



{A\m))M = {A^) + {5A\m))M. 



(C2) 



Hereafter we omit the subscript n for brevity. Tlie average {A^{m))M can be considered 
as a random variable with a variance We assume, that the measurement of A^{m) is 
a stationary ergodic process, and o'\j^^ — )■ 0. For a sufficiently large but finite number 
of diffraction patterns M satisfying ergodicity condition, the variance a\j can be estimated 
aPl 



M 



1 

M S 



(C3) 



Am=0 



Here C(Am) (Am = m — m') is the autocovariance of A'^{m) 



C{Am) = {A\m)A\m + Am)) - {Ay, 



(C4) 



where statistical averaging is performed over a large ensemble of diffraction patterns. 



Using the estimated value of crf^ [Eq. (C3)] we can determine a confidence interval for 
{A'^)m, for a given e. This can be done applying Tchebycheff 's inequality^. 



P{\{A')m-{A')\<6}>1 



a 



M 



(C5) 



Eq. (C5) determines the probability P{x} that the average {A'^)m lies in the interval {A'^)±e. 



As an example, we consider the system of clusters with a uniform distribution of 
orientations (pk and determine the confidence interval for {A'^)m, for a given accuracy e. 
Taking into account that for such a system the probability density function p[A^(m)] is 
knowrPl 

p[A''{m)] = ATe-^^'^'"), < A\m) < 1, (C6) 



we obtain for the autocovariance (C4) 



C(Am) 



Am = 0, 
0, Am ^ 0. 



Substituting Eq. (C7) into Eq. (C3) we have for the variance 

1 



a 
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(C7) 



(C8) 



Using this value of cr|^ we can determine the confidence interval for {A'^)m according to 
Eq. (C5). For the 10% accuracy e = 0.1(74^) = 0.1/A^, where we took into account that 



(A^) = l/N for a system considered here. With this we obtain from Tchebycheff inequality 



[Eq. (C5)] 



P (0.9(^2) < {A^)m < 1.1(^2)} > 1 



102 



(C9) 



From this inequality the probability P that {A'^)m lies in the interval from 0.9(^4^) to 1.1(^4^) 
is P > 0.9 for M = 10^ and P > 0.99 for M = 10^ diffraction patterns. 
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TABLE I. Fourier components of the CCFs (n = 10, 20) for a single centered pentagonal cluster 
l"'{q)] and for a completely oriented dense system of clusters (C^ '^'^'^^'^) described in Sec- 
(Fig. pi. The relative standard deviation is defined as 6C^ = (C" 



^(jn dilute 

tion 



IV 



dilute'^ 



n dilute^ /f^n dilute 



Model 



q, [nm- 



qi = 0.023 



(?2 = 0.029 



93 = 0.036 



(74 = 0.043 



55 = 0.059 



95 = 0.059 



dilute — 



£"(9), [Eq. (9)] 



0.225 



0.142 



0.028 



0.366 



0.143 



0.152 



c: 



'71 dense 



[Fig. 



d)] 



0.245 



0.101 



0.044 



0.313 



0.156 



0.175 



5Cl\ % 



-29 



57 



-14 



15 



TABLE IL Fourier components of the CCFs (n = 10,20) for a dilute ((C^ ^^^i"*^)) and dense 
(C" '^^^^'^'j systems composed of 121 centered pentagonal clusters with the Gaussian distribution of 



orientations described in Section 



IV 



(Fig. [4]). The parameter dC^ is defined as in Table |l| 



Model 




/^20 


q, [nm^"'"] 


91 = 0.023 


92 = 0.029 


93 = 0.036 


94 = 0.043 


95 = 0.059 


95 = 0.059 


dilute^ ^ [Eq 


■ 


121] 


0.152 


0.096 


0.019 


0.247 


0.096 


0.031 


(^n dense ^ [Fig. 


4 


d)] 


0.263 


0.078 


0.016 


0.233 


0.071 


0.024 


6C'g\ % 


73 


-18 


-15 


-6 


-26 


-23 
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TABLE III. Fourier components of the averaged CCFs {CJ^)a4 (n = 10,20) for a dilute 
((C^ '^''"*'^)m, M — ;> oo) and dense ((C^ '^^°'^^)m) systems composed of = 11 centered pen- 
tagonal clusters with the uniform distribution of orientations. The parameter 6Cg is defined as in 
Table HI 



Model 


(Cg°)M 




q, [nm^i] 


qi = 0.023 


(?2 = 0.029 


93 = 0.036 


94 = 0.043 


95 = 0.059 


95 = 0.059 


dilute^ [Eqs. ( 


221] 


2.05E-2 


1.29E-2 


0.25E-2 


3.33E-2 


1.30E-2 


1.38E-2 


^^ndensc^^^^ [Fig. 7 


a)] 


2.56E-2 


1.67E-2 


0.56E-2 


3.56E-2 


1.37E-2 


1.49E-2 


SC^, % 


25 


29 


120 


7 


15 


13 



TABLE IV. Fourier components of the CCFs (n = 10, 20) for the case of a fully coherent 
scattering [Cq'^°^] and a partially coherent scattering {Cq'^'^"^, ?coh = 600 nm) from a system of 
11 clusters with uniform distribution of orientations. The relative standard deviation is defined as 



SC" = [C" (7" cohj^^n. 



Model 






9, [nm^"'"] 


91 = 0.023 


92 = 0.029 


93 = 0.036 


94 = 0.043 


95 = 0.059 


95 = 0.059 


Cncoh^ [Fig. 6 


;c)] 


0.077 


0.050 


0.024 


0.107 


0.045 


0.022 


qipcoh^ [Fig. J 


?»] 


0.050 


0.040 


0.008 


0.063 


0.031 


0.019 




-35 


-20 


-68 


-42 


-31 


-12 
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Incident 
x-ray beam 
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Disordered 
2D sample 
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y 



y 



(a) 



(b) 



FIG. 1. (Color online) (a) Geometry of the diffraction experiment. A coherent x-ray beam illumi- 
nates a 2D disordered sample and produces a diffraction pattern on a detector. The direction of 
the incident beam is defined along the z axis of the coordinate system, (b) A disordered 2D sample 
composed of centered pentagonal clusters (an enlarged view of the cluster is shown in the inset). 
The case when all clusters have random position and orientation in the 2D plane is shown. 
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FIG. 2. (Color online) The structural terms |L"(g)p [see Eqs. ([7| and (A3)] calculated for a single 
centered pentagonal cluster as a function of q for (a) n = and (b) n = 10, 20 and 30. These 
terms are normalized by the value of the form- factor \f{q)\'^ of a sphere, (c) The normalized values 
C^{q) = \L^{q) / L^{q)\'^ of the same Fourier components. Arrows indicate the q values at which 
the CCFs presented in this paper were calculated (circles mark the values of the specific Fourier 
components at these positions). 
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FIG. 3. (Color online) Cross-correlation analysis of the diffraction data for a completely oriented 
dense system {{R)/d ~ 1.5) consisting of 121 pentagonal clusters, (a) Diffraction pattern; (b) 
Angular averaged intensity {I{q,ip))ip. (inset) Angular distribution of clusters in the sample. In 
this case all clusters have the same orientation; (c) CCFs Cq{A) calculated at the selected q values, 
indicated in (b); (d) Fourier spectra of the CCFs for each selected q value, (insets) Enlarged 
view of the first 25 Fourier components. 
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FIG. 4. (Color online) Cross-correlation analysis of the diffraction data (similar to Fig. [s]) for 
a system consisting of 121 pentagonal clusters described by the Gaussian distribution of cluster 
orientations with the standard deviation cr,^ = 0.05 • 360°/5 = 3.6° and zero mean [see inset in (b)]. 
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FIG. 5. (Color online) Cross-correlation analysis of the diffraction data (similar to Figs. [3] and [4]) 
for a system consisting of 121 pentagonal clusters described by the uniform distribution of cluster 
orientations in the angular range (— 7r/5,7r/5) [see inset in (b)]. 
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FIG. 6. (Color online) Diffraction patterns and Fourier components of the CCFs calculated for 
three disordered systems, containing 11 clusters in different, uniformly distributed orientations. 
Angular distribution of clusters in each sample is given in the angular diagrams. Fourier spectra 
were calculated for each system at three q values, q2 = 0.029 nm^^, qs = 0.036 nm~^, and 
94 = 0.043 nm-^ 
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FIG. 7. (Color online) Fourier spectra of the CCFs {Cg{A))M, averaged over M = 1000 diffraction 
patterns. Three disordered systems composed of (a) 11, (b) 60 and (c) 121 clusters with the uniform 
distribution of angular orientations are considered. Fourier spectra were calculated for each system 
at four q values, q2 = 0.029 nm~^, ^3 = 0.036 nm~^, = 0.043 nm^^ and q^ = 0.059 nm~^. Red 
solid line in each Fourier spectrum defines the magnitude of the inter-cluster contribution estimated 



using Eq. (BIO) 
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FIG. 8. (Color online) Diffraction patterns and Fourier components of the CCFs calculated for a 
disordered system consisting of 11 clusters in the same positions and orientations as in the sample 
3 in Fig. |6j Different values of the transverse coherence length /coh of the incoming x-ray beam: 
(a) ^coh = 1-2 /um, (b) /coh = 600 nm, and (c) /coh = 300 nm were considered. Fourier spectra were 
calculated for each case at three q values, q2 = 0.029 nm~^, q-^ = 0.036 nm~^, and ^4 = 0.043 nm~^. 
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(a) 11 clusters 



(b) 60 clusters 



(c) 121 cluster 
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FIG. 9. (Color online) (Color online) The evolution of different terms in the expansion of the Fourier 



component (C^ )m (n = 10) [Eq. (19)] as a function of the number M of diffraction patterns used 
in the averaging. The results are presented for the same q values as in Fig. [7j The following 
contributions are shown: (S'f )j\//(|/°(q')p)A/ (blue curve), (S'2 + S'^)M/{\I^{q)\'^)M (black curve), 
{S2) M / {\I^{q)\'^) M (red curve), and the sum of all terms {C^)m (green curve). Three disordered 
systems composed of (a) 11, (b) 60 and (c) 121 clusters with the uniform distribution of angular 
orientations are considered. 
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(a) (b) (c) 




FIG. 10. (Color online) Derivation of the RDF for a square sample, (a)-(c) Three cases when 
the reference particle is placed at the position Pi {i = 1,2,3), which is located from the sides of 
the square at distances u and u, are shown. Ppi^^t{R,u,i') is proportional to (a) the length of the 
circumference of the radius R, (b) the length of the single arc bounded by the square, and (c) the 
total length of three arcs bounded by the square (in all cases solid red parts of the circumference); 
(d) RDF for spherical particles of a diameter d = 600 nm distributed within the square 2D sample 
with the length of a side of the square / = 10 /um; Three model systems used in calculations of 
inter-particle contribution consist of (e) spherical particles (f) triangular clusters and (g) centered 
pentagonal clusters. 



44 




FIG. 11. (Color online) Fourier components of the average CCFs corresponding to the inter-particle 



correlations {C 



q (int—clust)' 



(for n = 4). Calculations were performed using asymptotic expressions 



[Eqs. (BIO), (B16a) and (B16b)] for a system composed of spherical particles (red circles). For 
comparison, the results of direct averaging of Eq. ( |B1| ) over the set of M = 1000 diffraction patterns 
for a system composed of spherical particles (blue diamonds), trigonal clusters (yellow triangles) 
and pentagonal clusters (black squares) are shown. The number of particles in the system (a) 
N = 11, (b) N = 60, and (c) N = 121. 
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